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ABSTRACT 



Context. The black hole at the Galactic Center, Sgr A*, is the prototype of a galactic nucleus at a very low level of activity. Its radio 
through submm-wave emission is known to come from a region close to the event horizon, however, the source of the emission is still 
under debate. A successful theory explaining the emission is based on a relativistic jet model scaled down from powerful quasars. 
Aims. We want to test the predictive power of this established jet model against newly available measurements of wavelength- 
dependent time lags and the size- wavelength structure in Sgr A*. 

Methods. Using all available closure amplitude VLBI data from different groups, we again derived the intrinsic wavelength-dependent 
size of Sgr A*. This allowed us to calculate the expected frequency-dependent time lags of radio flares, assuming a range of in- 
and outflow velocities. Moreover, we calculated the time lags expected in the previously published pressure-driven jet model. The 
predicted lags are then compared to radio monitoring observations at 22, 43, and 350 GHz. 

Results. The combination of time lags and size measurements imply a mildly relativistic outflow with bulk outflow speeds of yP ^ 
0.5 - 2. The newly measured time lags are reproduced well by the jet model without any major fine tuning. 

Conclusions. The results further strengthen the case for the cm-to-mm wave radio emission in Sgr A* as coming from a mildly 
relativistic jet-like outflow. The combination of radio time lag and VLBI closure amplitude measurements is a powerful new tool for 
assessing the flow speed and direction in Sgr A*. Future VLBI and time lag measurements over a range of wavelengths will reveal 
more information about Sgr A*, such as the existence of a jet nozzle, and measure the detailed velocity structure of a relativistic jet 
near its launching point for the first time. 
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1. Introduction 

■ The Galactic center hosts by far the best constrained super- 
' mas sive black hole c andidate: the compact radio source Sgr A* 
, (see iMelia & Falck3l200L for a review). Its mass is believed to 
' be aroi ind 4 X lO^Mp, based on stellar proper motion measure- 

■ ments (" Schodel et alj|2002l: iGhez et al.ll2005l) . Linear polariza- 
tion measurements indicate that it is ext reme ly underfed, with an 
accre t ion rate of less than 10"^Mn/vr (lAg ol 2000; Bower et al.l 
120051; iMacquart et al.ll2006t iMarrone et al.li2007.) . The accretion 
rate and low radio flux put Sgr A* at the tail end of the local 
luminosity function ("N agar et al.ll2005 ^ of low-luminosity active 
galactic nuclei (LLAGN). This makes Sgr A* an ideal laboratory 
to study supermassive black hole physics in the quasi-quiescent 
state in which most galactic nuclei exist today. 



Sgr A* has been det ected at radio ( Balick & Brown|[l974l) 
and no w near-infrared (Genzel et al.l 1200 3) and X-rav wave- 
lengths (Bagaiiofi' eLal. 2001). The radio spectrum of the source 
is variable, slightly inverted, and peaking i n a submm-bump 
which originates close to t he event h orizon ("Zvlka et al. 19921: 
falcke et al. 1998, 2000; M eUa & Fa lcke 2001; MiyazaMeFaD 
2004t lEckart et alj 1200 6). The latter is of particular impor- 



cast by the event horizon dPalcke et"ani2000l: iHuang et"al]|2007l ; 
iBroderick & Loebl '2006). However, until recently no structural 
information was available for Sgr A*. At wavelengths shorter 
than that of the submm-bump, the resolution of current tele- 
scopes is insufficient and at long wavelengths, where high- 
resolution very long baseline interferometry (VLBI) techniques 
can be used, the source structure is blurred by interstellar scat- 
tering. 

This ambiguity has led to a longstanding debate about the 
actual nature of the Sgr A* emission. One class of models sug- 
gests that the radio through X-ray emission is caused by ac- 
creting hot plasma flowing into the black hole (Meliai 119921 ; 
lNaravanetal]|1998l; lOuataert & Gruzinovll2000l) . On flie other 
hand, it has been suggested that Sgr A* resembles the com- 
pact radio cores of active galactic nuclei (AGN) and therefore 
most of the emis sion is associated with a (mildly) relativis- 
tic outflow or let (iRevnolds & McKeelll980l;lFalcke et al.lll993 



Falcke & BiermannI 119991; iFalcke & MarkoffI 120001; I Yuan et all 



200l) . 



tance since it may eventually allow imaging of the shadow 
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Only recently have measurements of the intrinsic size of 
Sgr A* become avail able (Bower etal. 2004; Shen et al. 2005]; 
iDoeleman et aLllToOSl) . providing crucial new input. The new in- 
trinsic size measurements agree well with the predictions of the 
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traditional jet model (iBower et alj 120041: iMarkoff et aT] l2007h . 
however, a direct confirmation of an outflow is still lacking. 

Clearly, additional information is required to determine the 
speed and direction of the flow responsible for the emission in 
Sgr A*. Such additional information has now become available 
with the first reliable time lag measurements of radio ou tbursts 
at diff'erent wavelengths ( Yusef- Zadeh et al.ll2006all2008l) .These 
observations show that high radio frequencies lead the lower ra- 
dio frequencies by some 20 minutes around 43 GHz. Because 
the radio emission is considered to be optically thick due to its 
flat-to-inverted spectrum, and the synchrotron loss timescale is 
much longer, the radio flux variations are tracing actual adiabatic 
expansion or contraction of the emitting plasma. This scenario 
is in marked contrast to observations in the optically thin part of 
the spectrum at near-infrared (NIR) and X-ray bands, where the 
cooling time scales are faster than adiabatic. At these higher fre- 
quen cies, observations (iMarrone et al.ll2008t iDodds-Eden et al.l 
I200&) show a near simultaneity between NIR and X-ray flares 
within minutes and a delay between X-ray/NIR with respect 
to the radio emission on ther order of hours. The expectation 
therefore is that radio timing observations trace bulk plasma 
properties, while X-ray/NIR variability is dominated by heat- 
ing and cooling of particle energy distributions in the plasma. 
Which physical parameters de termine a potential l ag between 
X-rays/NIR and radio/submm (iMarrone et al.ll2008l) . is not im- 
mediately obvious. 

In this paper we focus on the radio time lag data and size 
measurements to obtain information on the plasma flow speed. 
To do this we re-derive the intrinsic size of Sgr A* by combining 
all existing VLBI data in Sec. 12. 11 thereby resolving some of the 
apparent discrepancies between the results of different groups 
in the literature. We then compute the predicted time lags for 
various inflow/outflow speeds in Sec. 12.31 and present time lag 
predictions of the canonical jet model in Sec. 12.41 Here we also 
present the only analytical velocity profile of a pressure-driven 
jet in a closed form. The predictions are then compared with the 
data under the assumption that the region causing the variabil- 
ity roughly follows a similar size-frequency relation as seen by 
VLBI, tracing the bulk of the plasma. Our main conclusions are 
then summarized and discussed in Sec. [3] 

2. Size and time lag data in Sgr A* 

2.1. VLBI size of Sgr A* 

The radio size of Sgr A* is extremely difficult to determine 
for several reasons. The radio source itself is very compact and 
hence VLBI techniques have to be used, where radio telescopes 
with separations of several thousand kilometers are combined 
to obtain interferometric information of the source structure. 
However, the major high-frequency VLBI telescopes are in the 
Northern hemisphere, making Sgr A* a low-elevation source 
which is difficult to calibrate. Moreover, the source is located in 
the Galactic center behind a large scattering screen that broad- 
ens the intrinsic source size significantly at long wavelengths. 
To escape scattering eff'ects requires observing at shorter wave- 
lengths, which are even more difficult to calibrate . Hence, the 
breakthrough for the detection of the intrinsic size dBower et al.l 
|2004<) came via the in troduction of closure ampUtude analysis 
dPoeleman et al.ll2001h . a method relatively insensitive to com- 
mon station-based calibration errors. 

Closure amplitudes provide good means to measure the 
source size with very high accuracy, especially if the source 
structure is simple. Since the broadening of the source struc- 
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Fig. 1. Measured radio source size (major axis) of Sgr A* as 
function of observing wavelength in centimeters. 



ture by scattering follows a law (iDavies et al.l 
Ivan Langevelde et all [19921 iLo et all 1 19981 iBower etal] 
the actual source size <^sgi a* is given by 
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(1) 



where 0obs and 0scatt are the actually observed and the expected 
scattering size respectively, ^scatt can be obtained by measuring 
the source size at long wavelengths, where the intrinsic size is 
negligible, and extrapolating with a /l^-dependence to shorter 
wavelengths. The validity of thi s extrap olation of the A^-law has 
been discussed by iBower et al.l (l2004l) and demonstrated using 
the measured Gaussianity of the scattered image. 

In the following we employ this procedure using all currently 
available data, and discuss the origin of apparently conflicting 
results for the wavelength-size structure of Sgr A*. 

There are currently only four pap ers which contain reliable 
major and minor axes for Sgr A *: Bowe r et al.1 (|2004|) for A6 
cm to Al mm data. Bower et al.l (2006 ) for A2 4 cm to All cm 
data, and IShen et alj (12005) and iSheni (l2006h for A3 mm and 
Al mm data. Sizes at wavelengths longer than 20 cm are from 
VLBA closure-amplitude measurements and at longer wave- 
lengths high-quality VLA data is availabl e. There is also a clo- 
sure amplitude size at A3 mm from Poeleman et al J (12001), how- 
ever, that was only reliably obta ined for a c ircular Gaussian 
source and has been superseded bv lShenI (|2006|) . who fit an ellip- 
tical Gaussian. In addition there was one older measurement at 
/II. 3 mm, bv Krich baum et al.l d 19981) . based on a single baseline 
detection. Th e latter has now been su perceded by a more recent 
detection by iDoeleman et al.l ( l2008l) . based on three baselines 
and higher signal-to-noise ratio. The A1.3 mm observations yield 
the smallest sizes and the largest excursion from the scattering 
law. We therefore include these data points in our analysis for 
completeness, even though they do not represent a closure am- 
plitude measurement and one cannot distinguish between major 
and minor axis. Their inclusion, however, does not change our 
results significantly. 

Figure [1] shows the size data of Sgr A* as function of wave - 
length together with the scattering law from lBower et alj (l2006l) . 
0,scatt = (L31 + 0.02) mas {A/cm}^. One can clearly see how the 
overall size of Sgr A* follows the A^-law closely at long wave- 
lengths. 

The next step is to subtract the scattering law in quadrature 
from the observed size according to Equation [1] For this, the 
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exact normalization of the scattering law is vi tally important. 
The no rmalizations used by iBower et alj ( l2006l) and IShen et all 
( l2005h differ only slightly. This has little impact on the in- 
trinsic source size at A3 mm and A7 m m, but markedly af- 
fects the size at longer wavelengths. As iBower et al.l (l2006l) 
sho wed, this changes th e size-vs-wavelength relation (size oc 
A'"). IBower et all (|2006|) find power laws in the range between 
m - 1.3 and m — 1.7. dShen et a l. 2005), who use only short- 
wavelength data, find m = 1.09. 

The biggest problem, therefore, is the systematic uncertainty 
introduced by the inclusion or non-inclusion of long-wavelength 
data sets. We investigate this uncertainty in the following dis- 
cussion. Note however that the difference in the scattering law 
primarily affects intrinsic sizes at long wavelengths; short wave- 
length sizes are largely unaffected because the contribution of 
the scattering angle to the observed size is much less. 

2.2. Robustness of the Sgr A* size measurements 

Figure |2] shows the observed sizes divided by A^. Here we have 
averaged the data for the various observing bands, in order to 
avoid having the final fit be biased by the number of observations 
in one band. For the averaging we divided the sizes by to take 
out the frequency dependence, and weighted them by their error 
bars. This giv es us one data point per band. In particular, all 20 
cm data from iBower et al.l t2006.) are averaged into one point 
here. The error bars we show are the standard deviations of the 
measurements in one band, where multiple measurements were 
available. In principle this should be a more robust measure of 
the error. 

The non-homogeneous error distribution is problematic, but 
as it is a Umit of the available observational data base, it can- 
not be overcome. The three data points at A7 mm, A3. 5 cm, 
and A2Q cm tend to dominate any fitting and a combined multi- 
parameter fit of scattering-law and intrinsic size does not con- 
verge. Therefore it is customary to only fit the scattering law 
to long-wavelength data. The range of currently used scattering 
laws then depends exclusively on which data to include. Any un- 
known systematic error at A3. 5 cm or A20 cm would drastically 
affect the result. To quantify the robustness of the inferred sizes, 
we performed a series of weighted fits to the data below Al cm, 
with one random data point dropped. Doing this we find a range 
of possible scattering laws (Fig.|2]) given by 
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Fig. 2. Measured radio source size (major axis) of Sgr A* di- 
vided by A^ as function of observing wavelength in cm. The solid 
line represents the average scattering law used here, where the 
dashed lines indicated the 3cr limits found by randomly drop- 
ping one data point. 



we then fitted the intrinsic size law and determined the slope pa- 
rameter m. We find that the distribution of m is non-Gaussian, 
with a more extended tail towards smaller values. The median, 
however, is again at m - 1.44 - 0.19 + 0.16, where the errors 
are the 25% and 75% quantiles, respectively. These may be the 
more realistic error estimates than the ones from the simple ana- 
lytic fitting. 

To improve on this result in the future, more and better 
closure-amplitude size measurements need to be obtained at 
longer wavelengths, especially at A2 and A6 cm. 

In any case, the combined set of currently available data and 
the error analysis confirm previous conclusions that there is a 
wavelength-dependent photosphere in Sgr A* from a stratified 
medium. As expected for optically thick synchrotron radiation, 
the optical depth is indeed frequency dependent. This means 
that observations of Sgr A* at two different radio wavelengths 
provide information about two different spatial scales where the 
emission originates. 



(1.36 ±0.02) masx(.l/cm)^ 



(2) 



This i nclu des the bes t -fit sc attering laws used by iBower et alj 
( 120061) andHhenetaD ( l2005 l) within 3 cr limits, which have scal- 
ing factors of 1.31 + 0.02 and 1.39 ± 0.02 respectively. 

Subtraction of this scattering law in quadrature yields a 
slightly revised intrinsic size as shown in Figure [3] The sizes 
at A2 cm and A3. 5 cm are relatively sensitive to the scattering 
law and therefore "negative" source sizes are possible within 
the errors. Negative sizes are treated as lower limits around zero 
with the respective error bars. We fit the error-weighted intrinsic 
source size with a powerlaw function, yielding: 



*SgrA* - (0.52 + 0.03) mas x (A/cm) 



1.3±0.1 



(3) 



Again, this is consistent with the previous results and will 
be used in the following analysis. We have further verified this 
result by running a Monte Carlo simulation, excluding the 1.3 
mm data, by randomly varying the observed data and the scat- 
tering law within the errors quoted here. To each of these trials 



2.3. Variability and time lags 

In addition to the size measurement, we now have another new 
crucial parameter: the time lags between different wavelengths 
during small-scale variability outbursts. In the absence of direct 
imaging of source substructure, this provides the only means to 
determine flow or signal speeds in Sgr A*. 

The overall variability of Sgr A* has been established for a 
long time. The most comprehensive data sets stem from long- 
term monitor ing programs with the Green Bank In terferometer 
dFalckdl 19991) and the VLA ( IHerrnstein et alj|2004l) at cm wave- 
lengths. The reported rms variations of the radio spectrum are 
2.5%, 6%, 16%, 17%, and 2 1% at wavelen gths of 13, 3.6, 2, 
1.3, and 0.7 cm respectivelv. iMacquart & Bowed (12006) argue 
that most of the variation at longer timescales (several days) and 
at long wavelengths is due to interstellar scintillation. However, 
for time scales less than four days the variations may be intrinsic 
with an rms of ~ 10% for wavelengths 0.7-3 cm. V ariability is 
also seen at mm and sub-mm wavelengths dZvlka e t al. 1995; 
IZhao et al.1 l200l iMivazaki et al.1 12004 iMauerhan et al.. ,20051; 
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Fig. 3. Intrinsic radio source size (major axis) of Sgr A* obtained 
by subtracting the scattering diameter in quadrature. The soHd 
line represents the best fit powerlaw. The upper and lower dashed 
lines are the intrinsic sizes fitted by Bower et al. (2006) and Shen 
et al. (2005) respectively. 



iMarrone et al J2008|) with yet lar ger rms variations and outbursts 
of a factor of several over the quiescent level. 

In most cases, where multiple wavelengths were observed, 
the time coverage was not dense enough to find a reliable 
time lag between two wave lengths, despit e several attempts. 
Recen tly, lYusef-Zadeh et al.l (j2006a) and lYusef-Zadeh et alj 
(I2OO 8) published data obtained with the VLA in fast switching 
mode allowing quasi-simultaneous high-time resolution mea- 
surements of time variability in Sgr A* at two different wave- 
lengths. They find a lag between Al.3 and 0.7 cm on the or- 
der of 20 minutes . Takin g the weighted average of Table 1 in 
lYusef-Zadeh et al.l (l2008l) one finds a lag of 21 ±3 minutes. 
Millimeter and submm-millimeter wave timing observations by 
the same group are less significant, but seem to go in the same 
direction, with a lag between 22 and 350 GHz of 65 ± +10-23 
minutes ( Yusef-Zadeh et al. 2008). 

The sign of the lag between 43 and 22 GHz (^^0.7 & 1.3 
cm) is such that the shorter wavelengths lead the longer ones. 
Keeping in mind that shorter wavelength emission originates at 
smaller size regions, this immediately implies that bursts propa- 
gate outwards from small to larger scales. 

Given that we know the projected size s = ^SgrAt^^GC of 
Sgr A* from observations - Dqc is the Galactic center distance 
Z) = 8 kpc (Eisenhauer et al. 2003) - the time lag provides a 
straightforward estimate for the flow speed. Using equation Eq.[3] 
we find that the intrinsic size of Sgr A* is (pi = 0.73 mas and 
<p2 = 0.32 mas or sj = 8.8 x 10'^ cm and 52 = 3.9 x 10'^ cm 
at Al.3 cm and A7 mm respectively. Hence, As is ~ 27 light 
minutes. Given a time lag on the order of At = 20 min the flow 
velocity is v = (s\ - S2)/At - 1.4c. Therefore, Sgr A* needs to 
a harbor at least a mildly relativistic outflow, even if one allows 
for an error of ~50%. Projection effects would tend to increase 
this value even further 

Alternatively, if one has a model for a flow speed v(s) one 
can predict the time lags with Ar = (si/v(si) - S2/v(*2))- The 
easiest model is one with a constant flow speed v{s) =const. To 
allow for relativistic speeds we write this as v(s) - yfic, where 
y - - and yS - vjc. For this subsection, we will ignore 
projection effects for the sake of simplicity. The time lag then is 
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Fig. 4. Expected time lag of Sgr A* versus wavelength relative 
to /II. 3 cm (22 GHz) for the observed size-wavelength relation 
and a proper flow or signal speed of y/3 - 1 (red, solid line) 
or y/3 = 0.5 and 2 (orange, dashed). The data points are mea- 
sured time lags from Yusef-Zadeh et al. (2008). The top black 
solid line shows the Newtonian time lag for an outflow just at 
the escape speed. Long lags above that line would correspond to 
gravitationally bound outflows. 



At - Dcc^cpsgrA* /y/Sc. Figured shows the time lags for proper 
speeds y/3 in the range 0.5-2 c for the measured size-wavelength 
relation. 

We note that the source size in Sgr A* is close to linear with 
wavelength, hence, for constant velocity one expects a linear in- 
crease of the time lags with decreasing wavelength relative to 
a fixed reference wavelength. For comparison, we also show 
the time lags if the outflow would propagate always with the 
(Newtonian) escape speed ( y/2GM/r) for a 3.6 x 10^ Mq black 
hole. Here the time lags would become longer and grow non- 
linearly towards shorter wavelengths, since the escape speed is 
significantly slower than the speed of light. Figure|4]shows that a 
gravitationally bound flow would predict much longer time lags 
than what is actually observed. 



2.4. Time lags in tine jet model 

Given that the time lags suggest a mildly relativistic out- 
flow, it seems appropriate to investigate what a n actual jet 
model predicts. The basic jet model for Sgr A* (iFalcke et al.l 
1993': 'Falcke & Markoff '2000) naturaUy fits the spectrum, prop- 
erly predicted the low accretion rate of Sgr A* now inferred 
from polarization, a nd also was able to explain the VLBI size 
(iMarkoff^ et al.ll20d7 l). The only major property that could not be 
tested so far is in fact the flow speed. 

So far, the underlying assumption for the jet model has been 
that Sgr A* is not a strongly relativistic outflow. Energetically 
this is an optimal solution in terms of the ratio between tota l 
jet power and emitted synchrotron radiation ("Fa lcke et alJl993h . 
On the other hand, the sound speed for a relativistic plasma as 
well as the escape speed from the black hole are on the order of 
~ 0.5c, which sets a lower bound for a superson ic jet in Sgr A*. 

In the standard iBlandford & Konigll (1 19791) model for the 
flat-spectrum radio emission of compact jets, a constant veloc- 
ity is assumed and introduced as a free parameter. Falcke ( 199^ 
pointed out that this is in principle inconsistent, since the longi- 
tudinal pressure gradient would inevitably lead to some acceler- 
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ation of a modestly relativistic jet. As a first-order assumption 
the velocity was then assumed to be simply given by a purely 
pressure-driven wind in the supersonic regime. This approach 
had the advantage that the actual acceleration mechanism of the 
jet, which is likely magnetohydrodynamic in origin, could be 
treated as a black box. 

Simulations of the actual acceleration p rocess are actually 
very difficult and tim e consuming (e.g., iMeier e t al. 2001t 
iDe Vilhers et al.ll2005l) . However, they all start with some initial 
magnetohydrodynamic collimation regime (here referred to as 
the "nozzle"). After passing through the fast magnetosonic point, 
the flow eventually becomes over-pressured in a phase where the 
jet expands more or less freely into the ambient medium. The lat- 
ter situation is mainly add ressed by simulation s of pc-scale jets 
observed with VLB! (e.g..' Mimica et al.ll2008l) . 

Since for our simple Sgr A* jet model only the part af- 
ter the sonic point was considered, the only main parameter is 
then the location of the sonic point and the sound speed. For 
powerful, relativistic jets the sonic point is expected to be up 
to th ousands of Sc hwarzschild radii away from the central en- 
gine (iMarscher et a l. 2008) while for Sgr A* a relatively small 
value, of a few Schw arzschild radii, seems required by the data 
jMarkoff et al.ll20()7h . The magnetized plasma is here treated as 
a single-component fluid with adiabatic index 4/3 - i.e., in the 
relativistic limit of a "photon gas". The supersonic jet evolution 
is then calculated from the modified, relativistic Euler equation 
for a freely expanding jet propagating along the z axis in a cyUn- 
drical coordinate system, which we reproduce here from lFalckg 
(119961) : 



d I cj\ d 



(4) 



(A) - mpiic^ + Uj + Pj is the enthalpy density of the jet, Uj is 
the internal energy density, n is the particle density, and Pj - 
(F - l)Uj is the pressure in the jet (all in the local rest frame). 
With a "total equipartition" assumption one gets t/j ^ m-pnc^, 
hence co - il + T)U-^ and w/n = (1 + F)mpC^ =const at the sonic 
point z - Zo- In the free jet with conical shape the energy density 

evolves as J/j oc {y]fi]\ z"^ ■ The final equation is then given by 



Eq. 2 in'Falcke (1996). 

Note that for simplicity this equation lacks a gravitational 
term. This term becomes negligible quickly at larger distances 
and for unbound flows corresponding to typical radio frequen- 
cies for Sgr A*. This is clearly a deficiency when discussing the 
nozzle region in detail. In the following we will subsume this 
effect in the nozzle size as a free parameter. 

Previously the solution of the equation was only available 
numerically in the code. In in the following we present it in a 
closed form that allows testing it against time lag observations. 
For an adiabatic index of F = 4/3 the solution (see Fig. |5]) is 
given implicitly as 



8 + 121^) -281n(— )+421n(z) 



(5) 



where f'iy) = x is the inverse function of / such that f{x) - y, 
X - yjySj, and 



f(x) = -28 In jc -I- -42^/3x5/3 + 42^2. 



(6) 




7j and ySjc are the bulk jet Lorentz factor and velocity, z - Zjznozz 
is the dimensionless length along the jet axis (Z), and Znozz marks 



log z/Znozz 

Fig. 5. The proper flow speed (Eq. [S]) of a pressure driven jet 
plotted versus the logarithm of the distance, in units of the nozzle 
size Znozz, along the z-axis. 



the location of the jet nozzle. The equation thus has a critical 
point at z = 1 , where yj^Sj equals the proper sound speed y^fi^ - 

-\/r(r - 1)(F + 1)"', and is only valid in the supersonic regime 
z> 1. 

This relatively simple quasi-analytical description had first 
been developed for M81* and was then integrated into the Sgr 
A* jet papers thereafter. While naturally overly simplified, we 
retain it here, treating it as a published prediction. However, one 
should not consider this as the only possible solution, but rather 
as representative of a broader class of models for modestly rela- 
tivistic accelerating jets. 

Using this description we now calculate the time lags based 
on the assumption that any flare is essentially due to an increase 
in the accretion power This increased accretion will turn into 
an increased outflow rate, based on the "jet-disk symbiosis"- 
ansatz of a linear coupl ing between inflow and outflow rate 
dFalcke & Biermannlll995l) . The increased power and mass flow 
will then propagate along the jet essentially with the local flow 
speed. Here we ignore the slightly increased acceleration due to 
the increased longitudinal pressure gradient in an overdense re- 
gion, which would be a second order effect. 

We note that earlier we have argued that the X-ray flares in 
Sgr A* are not due to a similar increase in accretion, but rather 
due to additiona l heating or acceler ation of the internal parti- 
cle distributions (jMarkoff et al.ll200"l]) . This is entirely consistent 
with our approach here, since in the same paper we showed that 
such heating processes only marginally affect the radio flux in 
the optically thick region. Hence simultaneous radio-X-ray flare 
are not necessarily required. Radio flares, however, required an 
actual increase in accretion rate as also argued here. Of course, it 
is not inconceivable that a sudden increase in accretion rate also 
leads to additional heating and particle accreleration in the inner 
region of disk and jet. 

Indeed, recent observations iMarrone et al.l (l2008l) seem to 
show that there is a rather long lag (on the order of hours) be- 
tween X-ray/IR-flares and radio flares. This time scale is much 
longer than free-fall or rotational time scales and consistent with 
viscous processes in the accretion flow linking the two types of 
flares. 
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Fig. 6. Expected time lag of Sgr A* versus wavelength relative 
to /II .3 cm (22 GHz) for the observed size-wavelength relation 
and a flow speed according to the jet model for three inclination 
angles. The data points are the same as in Fig.|4] 



The predicted radio time lags in the jet model are then cal- 
culated as 



At = 



A(/>Dgc 
sin iPj{z)c 



(l -ySj cos . 



A0 = <pSgiA*(A)) - 0SgrA*(/l) and Ao - 1.35 cm is chosen as the 
reference wavelength. This formulation recovers the well-known 
formula for apparent superluminal motion (jSapp = A0D/Ar), if 
the implied flares were observed as moving blobs. 

For the dimensionless length we take z - 
<^SgrA*(/t)/0SgrA.(/inozz) with A^^zz = 0.8mm. The latter 
represents the next observing band above the highest currently 
available VLBI measurements and corresponds to a projected 
size of about 4Rg (Rg = GM,/c^), for a black hole mass of 
M, = 3.6 X 10^ M(7). This is also the typical nozzle size used i n 
specti-al fits (e.g., iFalcke & Markoffll200(]t iMarkoff et al.ll2007h . 

Figure |6] shows the expected time lag for the measured size 
and the velocity field of the pressure-driven jet. The prediction 
is consistent with the 21 min time lag found between A7 mm and 
A1.3 cm. We stress that this is based solely on the combination 
of the observed sizes and the previously published velocity field 
for the jet model. 

Quite noticeable is the quick rise of the time lag, relative 
to Al.3 cm, towards shorter wavelengths. The rise comes from 
the fact that the jet first needs to accelerate beyond the noz- 
zle, which yields initially slower flow velocities and accordingly 
longer time lags. This should be a characteristic signal of a noz- 
zle, which future time lag measurements could help to identify. 
Of course, one has to bear in mind that the model is overly sim- 
plistic and in reality this feature may look less drastic. In par- 
ticular, general relativistic effects will start to play an important 
role. Also, the exact location of this kink is very sensitive to the 
size of the nozzle, which is a free parameter within a factor of 
two or so and therefore difficult to predict. On the other hand, 
an exact localization of the kink would eff'ectively constrain the 
nozzle size. 

For future use, we also extrapolate the predicted time lag to 
longer wavelengths (Fig. [T). One can see that the lag becomes 
on the order of a day at cm wavelengths. This may therefore be 
difficult to observe, given the limited observability of Sgr A* in 
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Fig. 7. Same as Fig. |6] but extended to longer wavelengths. For 
reference we also show the time lags for a marginally bound 
outflow as in Fig.|4] 
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Fig. 8. Same as Fig. |7]but now we show the time lags for a 
jet model inclined by 60° for the different size-wavelength laws 
considered. This has a significant effect at long cm waves. 

the Northern hemisphere, but would certainly provide crucial in- 
formation on the large-scale structure of Sgr A* that is otherwise 
impossible to obtain due to the strong scatter broadening. 

In addition we consider the effect of the range of possible 
size-wavelength relations for Sgr A* in Figure [8] Not surpris- 
ingly the time lags do not show much of a difference at short 
wavelengths, but differ markedly at long wavelengths. 

3. Summary and discussion 

By now Sgr A* is probably the best studied supermassive black 
hole with imaging and timing information on scales very close to 
the event horizon over a wide range of frequencies. New VLBI 
measurements, which we have here revisited, have confirmed 
theoretical predictions that Sgr A* has a frequency-dependent 
photosphere at radio wavelengths, with sizes scaling roughly as 

^1.3±0.1^ 

The time lag of individual bursts seen at different wave- 
lengths provides a powerful new tool to constrain the physics 
at work in Sgr A*. Combining this timing data with the increas- 
ingly better intrinsic size measurements obtained with VLBI can 
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significantly constrain flow speeds. The latter is otherwise not 
measurable with direct imaging due to the extreme scatter broad- 
ening in the Galactic center 

The well-esta blished lag of ~ 20 min between Almm and 
/II. 3 cm found by lYusef-Zadeh et alj (l2006b) . together with the 
intrinsic size difference of ~ 27 light minutes at these wave- 
lengths, already suggests that the radio emitting plasma is un- 
bound and flows out with mildly relativistic s peeds. 

Thi s is in contrast with the conclusions bv lYusef-Zadeh et"an 
(1200 6b') who derive a subsonic and sub-relativistic outflow from 
the same time lag data. However, the authors base t heir conclu- 
sion s solely on a fit of their light curve to a simple van der LaanI 
(11966.) model without ever considering the VLBI size measure- 
ments. As mentioned by these authors, the van der Laan model 
describes the adiabatic expansion of a spherical plasma blob and 
cannot describe bulk outflow, which we now know is a major fac- 
tor in the extragalactic radio sources it was developed for. The 
van der Laan mod el predicts source sizes of 8 Schwarzschild 
radii (e.g.. Fig. 3 in lYusef-Zadeh et alj|2006bl) at 22 GHz while 
the measured radio size at 22 GHz is 80 Schwarzschild radii. 
The variations of ^ 20% of the total flux would thus have to 
be produced by a region that is only ~0. 1 % of the total volume 
compared to the bulk of the plasma. We find this scenario un- 
likely. 

The pressure-drivenjet model dFalcke & Markofil2000h . that 
has been successfully used to fit size and spectrum of Sgr A* al- 
ready, is quite consistent with the combined size and time lag 
data. No particular adjusting of parameters is necessary with re- 
spect to published jet models. The main free parameters are the 
nozzle size and the inclination angle, for which we have used 
canonical values. Coincidentally, both parameters do not affect 
the A7 mm-/il.3 cm lag very much. However, as one can see in 
the figures, these parameters will become important once time 
lags at other wavelengths are available. 

The sensitivity to parameters can be turned around to state 
that measurements at other wavelengths in the future will pro- 
vide invaluable information about the structure of Sgr A*. If time 
lags can be found at A2 cm or A3.5 cm, this would constrain in- 
clination angle and the size-wavelength relation much better and 
even more clearly distinguish between models. 

Moreover, time lags at A3 mm and Al mm could start to 
show evidence for the acceleration region of the outflow ("the 
nozzle"), which would be a unique diagnostic for jet and accre- 
tion physics. However, here one has to caution that our simple 
analytic treatment naturally breaks down near the nozzle region 
and in the vicinity of the black hole. More so phisticated numer- 
ical magnetohydrodynamic calculations (e.g., iHawlev & Krohkl 
12006 ) are clearly needed to understand the submm-wave emis- 
sion. Also, since the jet likely originates from an inflow some- 
where close to the nozzle region, this inflo w could also con - 
tribute to the submm-bump emission (e.g.. ' Yuan etai]|2002h . 
These additional emission components might eventually de- 
crease the coherence of outbursts across wavelength in the 
submm & mm-wave region. 

So far we only have a single reliable time lag from Al mm 
to /II. 3 cm on which to base conclusions, with corroborating ev- 
idence from a tentative A 0.8 mm to A 1.3 cm lag. Therefore fur- 
ther time lag measurements at these and other wavelengths and 
VLBI closure amplitude measurements at long wavelengths are 
highly desirable. The former is challenging since the expected 
time scales at longer wavelengths are on the order of a typical 
observing run. The latter is challenging because typical VLBI ar- 
rays resolve out Sgr A* at long wavelengths and long baselines. 
Nonetheless, these observations are certainly worth the effort. 



They promise a wealth of detailed information about how plasma 
behaves very close to the event horizon of a supermassive black 
hole. Such coordinated multi-wavelength studies will eventually 
allow a range of more sophisticated models to be tested in great 
detail. 
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